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C/3 1 Abstract 



We study the Hawking effect in terms of the geometric phase acquired by a two-level atom 
as a result of coupling to vacuum fluctuations outside a Schwarzschild black hole in a gedanken 
experiment. We treat the atom in interaction with a bath of fluctuating quantized massless scalar 
fields as an open quantum system, whose dynamics is governed by a master equation obtained by 
tracing over the field degrees of freedom. The nonunitary effects of this system are examined by 
analyzing the geometric phase for the Boulware, Unruh and Hartle-Hawking vacua respectively. 
We find, for all the three cases, that the geometric phase of the atom turns out to be affected 
by the space-time curvature which backscatters the vacuum field modes. In both the Unruh and 
Hartle-Hawking vacua, the geometric phase exhibits similar behaviors as if there were thermal 
radiation at the Hawking temperature from the black hole. So, a measurement of the change of 
the geometric phase as opposed to that in a flat space-time can in principle reveal the existence of 
the Hawking radiation. 
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I. INTRODUCTION 



In quantum sense, a black hole is not completely black, but emits thermal radiation 
with a black body spectrum. This is the well-known Hawking radiation which is first found 
by Hawking in the paradigm of quantum field theory in curved space-time [1]. Ever since 
Hawking's original work, this striking effect has attracted a great deal of interest, and 
extensive work has been done trying to re-derive and understand it in various different 
contexts [1-10], ranging from Euclidean quantum gravity [2] to quantum entanglement [10]. 

In this paper, we will try to understand, in the framework of open quantum systems, the 
Hawking effect by studying the geometric phase acquired by a static two-level atom due to 
its coupling to the vacuum fluctuations outside a Schwarzschild black hole in a gedanken 
experiment. The geometric phase is an important concept in quantum theory. In 1984, 
Berry showed that when the Hamiltonian of a closed quantum system evolves adiabatically 
in a cyclic way, the state of the system acquires an additional phase besides the familiar 
dynamical one, which has a purely geometric origin [11]. Berry's work was soon generalized 
to nonadiabatic [12] and noncyclic evolution [13]. The geometric phase has been extensively 
studied both theoretically and experimentally [1 1]. In recent years, it has been demonstrated 
that geometric phase may have many potential applications, for instance, fault-tolerant 
quantum computation [15]. However, due to the inevitable interactions between the qubits 
and the environment, a pure state will be driven to a mixed one. Therefore, the study 
of geometric phase of systems under nonunitary evolutions becomes important, and some 
work has been done to generalize the geometric phase to open quantum systems [16-20], 
and the effects of different kinds of decoherence sources on the geometric phase have been 
analyzed [21-25]. By studying the geometric phase of an open system, we can draw some 
information about the nonunitary evolution of the system and the characteristics of the 
environment the system coupled to. In principle, the geometric phase factor can be measured 
by interfering the system which has undergone the above evolution with one that does not. 
Noteworthily, a new type of experiment on the quantum geometric phase of an open system 
undergoing nonunitary evolution has recently been done, where the phase is determined by 
measuring the decoherence factor of the off-diagonal elements of the reduced density matrix 
of the system [26] . 

Recently, Martin-Martinez et al. [27] have studied the geometric phase of an accelerated 
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atom which couples to a single-mode of a scalar field in vacuum and found that the phase 
differs from that of an inertial one as a direct consequence of the Unruh effect. Thus, in prin- 
ciple, an interferometry experiment can be performed to reveal the effects of acceleration. 
Later, we have considered a more realistic case, i.e., an accelerated two- level system which 
couples to all vacuum modes of electromagnetic (rather than scalar ) fields in a realistic mul- 
tipolar coupling scheme, and proposed using the geometric phase of non-unitary evolution 
to detect the Unruh effect [28]. In the current paper, we study the geometric phase acquired 
by a two-level atom interacting with a bath of fluctuating quantized massless scalar fields in 
vacuum outside a Schwarzschild black hole in a gedanken experiment, in the hope of gaining 
understanding of the Hawking effect in a new perspective. In the gedanken experiment, we 
first prepare the atom in a pure state, adiabatically put it at a fixed radial distance outside 
the black hole, then the atom interacts with the fluctuating scalar fields in vacuum, and as 
a result, it evolves from a pure state to a mixed state. During this nonunitary evolution, a 
geometric phase is acquired. The atom in the gedanken experiment can be thought of as 
being held static by an imaginary string with tension, or the the buoyancy of the Hawking 
radiation, or even some other unknown mechanism. The reduced dynamics of the atom is 
studied by tracing over the field degrees of freedom from the total system. The Boulware, 
Unruh and Hartle-Hawking vacua will be considered. At this point, let us note that the 
theory of open quantum system has been fruitfully applied to understand the Unruh [29], 
Hawking [9], and Gibbons-Hawking [30] effects respectively in a way that is different from 
the traditional one. 



II. THE MASTER EQUATION 

The Hamiltonian of the whole system takes the form 

H = H s + H^ + H' . (1) 
Here H s is the Hamiltonian of the atom, and its general form is taken to be 



i=i 

in which crj is the Pauli matrix, and uq is the energy level spacing of the atom. For simplicity, 
we can write it as H s = ^uoa-s. H<p is the Hamiltonian of the scalar field, whose explicit 
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expression is not required here. We couple the atom and the scalar field by analogy with 
the electric dipole interaction [31] 



#' = //<7 2 $(i,x) , (3) 

where \x is the coupling constant, which we assume to be small. 

The initial state of the whole system is characterized by the total density matrix p tot = 
p(0) <S> |0)(0|, in which p(0) is the initial reduced density matrix of the atom, and |0) is the 
vacuum state of the field. In the frame of the atom, the evolution in the proper time r of 
the total density matrix p tot satisfies 

^ 1 = -- h [H, Ptot (r)}. (4) 

The evolution of the reduced density matrix p(r) can, in the limit of weak coupling, be 
written in the Kossakowski-Lindblad form [29, 32-34] 

= ~[H eS , p(r)] + £\p(r)} , (5) 

where 

1 3 

The matrix and the effective Hamiltonian H c g are determined by the Fourier and Hilbert 
transforms of the correlation functions 

G + (x-y) = p 2 (0mxMy)\0) , (7) 

which are defined as follows 

POO 

0(A) = / dre lXT G + (r) , (8) 



oo 



KW= rr*m, (9) 

in which P denotes principal value. Then the coefficients of the Kossakowski matrix can 
be expressed as [29] 

dij = A 5ij - iB e ijk n k + C rurij , (10) 

where 

A=i[0(w o ) + 0(-w o )] , 5 = ^„)-M], C = -A. (11) 
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The effective Hamiltonian H c r contains a correction term, the so-called Lamb shift, and one 
can show that it can be obtained by replacing uq in H s with a renormalized energy level 
spacing Q as follows [29] 

H cS = ~hQa 3 = ^{u + l - [/C(-wo) - K{uo)\} <x 3 ■ (12) 

For convenience, the density matrix p is expressed in terms of the Pauli matrices, 



^ 1 = 1 



PiWtn . (13) 



Taking Eq. (13) into Eq. (5), one can show that the components of the Bloch vector |p(r)) 
satisfy 

§p\p{T))=-2H\p{T)) + \ri), (14) 
where \r)) = {0, 0, — 4B}, and Ti takes the form 

( 2A + C fi/2 \ 



(15) 



-ft/2 2A + C 
y 2iy 

Now we assume that the initial state of the atom is |^>(0)) = cos ||+) + sin ||— ), and then it 
evolves from a pure state to a mixed state due to interaction with the external environment, 
namely, a bath of fluctuating vacuum scalar fields. The time- dependent reduced density 
matrix of the atom can be easily worked out as 

p(r)= I 2 + 2A 1 j 2 I . (16) 

l e -2(2A+C)r+inr gin J _ g -4Ar cog 2 | _ ^(e" 4 ^ - 1) 

III. GEOMETRIC PHASE IN OPEN TWO-LEVEL SYSTEM 

In this section, we use the master equation governing the reduced dynamics of the atom 
in the last section to calculate the geometric phase of a static atom in the background of a 
Schwarzschild space-time, which is described by the following metric 

ds> ={l- 2 ~f) - T -^ : - r\d9* + sin 2 9d^) . (17) 

When a curved space-time in concerned, a delicate issue is how to determine the vacuum 
state of the quantum fields. Normally, a vacuum state is associated with non-occupation 



of positive frequency modes. However, the positive frequency of field modes is defined 
with respect to the time coordinate. Therefore, to define positive frequency, one has to 
first specify a definition of time. In a spherically symmetric black hole background, three 
different vacuum states, i.e., the Boulware [35], Unruh [36], and Hartle-Hawking [37] vacuum 
states, have been defined, each corresponding to a different choice of time coordinate. In 
the following, we examine the geometric phase of the atoms in all these vacua. 

A. The Boulware vacuum 

Let us begin our discussion with the Boulware vacuum, which is defined by requiring 
normal modes to be positive frequency with respect to the Killing vector d/dt. In order to 
analyze the geometric phase of the system, we make use of the master equation Eq. (5) and 
see how the phase is related to it. 

The geometric phase for a mixed state under nonunitary evolution is defined as [19] 



\fc=i / 
where Afe(r) and |</>a;(t)) are the eigenvalues and eigenvectors of the reduced density matrix 
p(r). In order to get the geometric phase, we first calculate the eigenvalues of the density 
matrix (16) to find A±(r) = \{l ± 77) , in which r] = \/ p\ + e -^ 2A+c ^> T sin 2 6 and p 3 = 
e -4Ar cos g _|_ B(^ e -4AT _ j£ j s obvious that A_(0) = 0. As a result, the only contribution 
comes from the eigenvector corresponding to A + 




(18) 




(19) 



where 




(20) 



The geometric phase can be calculated directly using Eq. (18) 




(21) 



Plugging the explicit form of 8 T , rj and into Eq. (21), we have 



I 



T 




R- Re AAT + cos0 



) 



1 



sin 2 fl + (.R-.Re^ + cos^) 2 



Vtdr . 



(22) 
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The result of this integral is rather tedious. Here, for small 70/wo) we perform a series 
expansion of the geometric phase for a single quasi-cycle and obtain to the first order [28] 

2vr 2 



7 « — 7r(l — COS 1 



(JO 



-(2B-Ccos6) sin 2 6 . 



(23) 



It is obvious that the geometric phase of the system is characterized by the coefficients B 
and C which determine the Kossakowski matrix a^-, and we now calculate them explicitly. 

Let us now calculate the geometric phase of a two-level atom which is held static at a 
fixed radial distance by some additional force, an imaginary string with tension, for instance, 
outside a Schwarzschild black hole. The trajectory can be described by 

1 



t(r) 



(t-t ), r(r)=r , 



(t) = 9n 



(r) = 0o • 



(24) 



In the Boulware vacuum, the Wightman function for the scalar field on the trajectory 
Eq. (24) is given by [38-40] 

roo —iujAt 1 , 

1 Y lm (0, (J)) I 2 [| Ri(u, r) | 2 + I Ri(u, r) | 2 ] du . (25) 



G + (x, x') 



lm 



lm \ 



The Fourier transform with respect to the proper time is 

] ^rf [I i^Vd^o, r) | 2 + I %(\^- , r) \ 2 ]9(X) , (26) 



° e iXr G + (x, x')dr = fi 2 ^ 2 

■DO j 



1tt\ 



where 9 (A) is the step function, and we have used the relation 

/ 



m=—l 



21 + 1 
An 



(27) 



The above Fourier transform is hard to evaluate, since we do not know the exact form of the 
radial functions Ri(oj, r) and R~i(u, r). Here, we choose to compute it both close to the event 
horizon and at infinity, with the help of the properties of the radial functions in asymptotic 
regions [40]: 

r 4w 2 

7TW > r 2M , 

^(2/ + l)|^(w,r)| 2 ~ \ i J (28) 



1=0 



1^(2/ + 1)| ^ M| 2 , r->oo 
1=0 

( 1 00 

Y> + i)| Bl („m 



4M 2 
4w 2 



1 _ 2M ' 

^ 1 r 



r — )■ 00 . 



(29) 



Plugging Eq. (26) into Eq (11), and using Eq. (28) and (29), we have 

A = B=-C^^[1 + f(u ^, r)] , (30) 
in which f(u,r) is defined as 

tt^r) = 7 — 2 J2W + V\Bdu)\\ (31) 
z=o 

and 70 = /i 2 wo/27r is the spontaneous emission rate. So, after some straightforward calcula- 
tions, we obtain, in the asymptotic regions, 

{ cos 9) + tt 2 -^ sin 2 9(2 + cos 9) (l + g 00 f^oy/gw, 2M)) r 2M , 

7 - 2 ; o ° (32) 

- 7r(l - cos 6 1 ) + it 2 - — sin 0(2 + cos 0) (l + 5(00 /(^oa/^oo, r)) r ->■ 00 . 

The result above shows that, after a single quasi-cycle, the static atom acquires a geo- 
metric phase, which depends on the geometric properties of its state space. The first term 
in the above equation is the one we would have obtained if the system were isolated from 
the environment, and the second term is the correction induced by the interaction between 
the atom and the environment. It is interesting to note that, in both asymptotic regions, 
the correction of the geometric phase takes the same form, which is what one expects in the 
Minkowski vacuum[24, 28] multiplied by a factor (1 + goo f(^oy/goo, r)). The additional cor- 
rection as opposed to the Minkowski vacuum, which is proportional to goo f(u>oy/g~oo, r ), can 
be attributed to the backscattering of the vacuum field modes off the space-time curvature 
of the black hole. This is similar in a sense to the case of the reflection of the field modes 
at a reflecting boundary in a flat space-time. From the definition of the grey-body factor 
f(oj,r) (31) we can see that, in the limit of r — > 00, g 00 /(^o\/#oo5 r ) tends to be zero, and 
the result returns to the case in the Minkowski vacuum, which means that the Boulware 
vacuum corresponds to our familiar concept of a vacuum state at large radii. 

B. The Unruh vacuum 

Now let us go on our discussion to the Unruh vacuum, which is supposed to be the 
vacuum state best approximating the state following the gravitational collapse of a massive 
body. The Wightman function in the Unruh vacuum is given by [38-40] 



/oo -iuiAt r I /, , „\ 1 2 , 

-00 L 



duj , (33) 



where k = 1/4M is the surface gravity of the black hole. Its Fourier transform is 



0(A) 



e lXr G + (x,x')dr 



SttX 



oo 

2 00 



Z=0 L 



6(Xy/g^)(l + 2l)\%(\ 



-4 



5-00, 



2Z)|ifr(AVw) 



(34) 



in which K r is defined as ft/v^oo • Inserting Eq. (34) into Eq. (11), one finds that 



and 



A = -C*^[l + g 00 /(c^Vflto, 2M) + T ] , 

7o 



-fi , oo/(^ov / ^oO' r )] 



(35) 



7o 



-4 = -C « — [1 + g 00 f(uj 0y /g^, r) + 



r — > 00 : 



7o 



B ~ — [1 + fi-oo f(uJ y/g^, r)] . 



(36) 



Straightforward calculations lead to 



tt(1 -cos#) + tt 2 — sin 2 0x 
2w 



7 



(2 + cos 0) (1 + 3oo f(oooV9^, 2M)) 
.2 7o 



cos 9 



(37) 



7T( 1 — COS 1 



7T 



2Wn 



sin x 



(2 + cos 9) (l + 5-00 /(woa/^oo^ r )) + #00 fi^oVOm, r 



cos 



00 



Form the equations above, we can see that, in the vicinity of the event horizon, the first 
term in the environment induced geometric phase is the same as the one in the Boulware 
vacuum, and the second one is proportional to a Planckian factor. At infinity, the Planckian 
factor is modified by a grey-body factor 300 f{^oy/9oo, r ), which is caused by the backscat- 
tering off the space-time curvature. In order to have a better understanding of the features 
of the geometric phase and of what it tells us about, it would be rewarding to examine 
here the geometric phase of a static atom in a thermal bath at temperature 1/(3 in a flat 
space-time. In this case, the Wightman function reads 



G + (x,x') 



— y 



m=— 00 



(t — f — im/3 - is)' 



and similar calculations leads to 
7 ~ — 7r(l — cos 9) 



7r 2 ^sin 2 0( 2 + cos0 



2u 



r'\2 



cos 9 



(38) 



(39) 



Comparing Eq. (37) and Eq. (39), one finds that, apart from the backscattering of the space- 
time curvature, the result is the same as that in the case of a two-atom system immersed 
in a thermal bath at an effective temperature T = K r /2n = T H /^fg^ in a flat space-time, 
where T# = k/2% is the Hawking temperature [1]. This clearly suggests that at the horizon, 
there is a thermal flux going outwards. However, when r — > 2M, the temperature T is 
divergent, since this effective temperature is a joint effect of both the thermal flux from 
the black hole and the Unruh effect due to that the system is accelerating with respect 
to the local free-falling inertial frame so as to maintain at a fixed distance from the black 
hole, and the acceleration diverges at the horizon. As the atoms are placed farther away 
from the black hole, the thermal radiation becomes weaker due to the backscattering off the 
space-time curvature. At the infinity, the grey-body factor vanishes and it returns to the 
Minkowski vacuum case. This suggests that, in the Unruh vacuum, no thermal radiation 
is felt at the infinity due to the backscattering of the outgoing thermal radiation off the 
space-time curvature. 



C. The Hartle-Hawking vacuum 



Finally, let us come to the Hartle-Hawking vacuum. The Wightman function in this 
vacuum is given by [38-40] 



Y, 



lm\ 



Aixu 



„—iu)At i AuiAt 

l Ri(uj,r) 2 1 



J g — 2ttuj/k 



2ttui/k 1 



du 



and its Fourier transform is 
0(A) 



/oo ^ 
e lXT G + (x jX ')dT = /2 2 J2 
■°° 1=0 



;i + 2Q 

8ttA 



%{\^r)\ 2 , \%{-\ 



9oo, r) 



1 2 n 



^ g — 2ttX/ K r 



^ g — 2ttX/ K r 



Similar calculations then lead to 



r^2M: 



_ 2TTUJ /n r _|_ 1 

5«^[l + <7oo f(u 0y /gVo,2M)\ , 



(40) 



(41) 



r — > oo 



A = -C 



To 



To 



e 2iru} /K r _|_ 



4 g27roJo/K r 1 



[1 + #00 f(^oy/goo, r)] , 



5 « — [1 + g 00 /(wov^oo. r )] 



(42) 



10 



and 



7 



-tt(1 -cos#) + tt 2 — sin 2 #x 

(2 + cos 9) (1 + ^oo /(woVflto, 2M )) + 



T (l + /(w 0v ^,2M))co S e 



?rfl 



COS I 



7T 



„2 7o 



2w n 



sin # x 



(43) 



(2 + cos 9) (l + 5( o /(w OA /^oo, r)) + 



-(1 + /(w 0v / ^00' r )) cos6) 



r — > oo 



Here, unlike that in the Unruh vacuum case, the geometric phase takes the same form in 
two asymptotic regions. The first term is the same as that in the Boulware vacuum, and the 
second term is composed by two parts, one is proportional to the standard Planckian factor 
and the other is proportional to a Planckian factor modified by a grey-body factor caused 
by the backscattering off the space-time curvature. This suggests that there are thermal 
radiation outgoing from the horizon and that incoming from infinity, both of which are 
weakened by the backscattering off the curvature on their way. In both asymptotic regions, 
the impact of the black hole to the geometric phase is the same as that of a thermal bath 
at temperature T = Tn/^g^ in a flat space-time. As r — > oo, the effective temperature 
becomes the Hawking temperature, since the acceleration needed to maintain the two-level 
atom system at a fixed radial distance vanishes, and the temperature is purely due to the 
thermal bath the black hole is immersed in. Therefore, the Hartle-Hawking vacuum does 
not correspond to our familiar concept of vacuum. It is actually a state that describes a 
black hole in equilibrium with an infinite sea of blackbody radiation [40]. 



IV. CONCLUSION 



We have studied the evolution of a two-level atom interacting with a bath of fluctuating 
quantum scalar fields at a fixed radial distance outside a Schwarzschild black hole in the 
framework of open quantum systems, and calculated the geometric phase of the atom in the 
Boulware, Unruh and Hartle-Hawking vacua respectively. We find, for all the three cases, 
that the geometric phase of the atom turns out to be affected by space-time curvature which 
backscatters the vacuum field modes. In the Unruh vacuum, the geometric phase behaves as 
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if there were an outgoing flux of thermal radiation which is backscattered by the space-time 
curvature. In the Hartle-Hawking vacuum, the geometric phase exhibits behaviors as if the 
atom were in a thermal bath at an effective temperature which reduces to the Hawking 
temperature at infinity. 

Generally speaking, the geometric phase is of significance only if we compare the situation 
of the atom to some other situations. For our current case, it is natural to compare the atom 
kept at a fixed position to the one released free. However, the calculation of the geometric 
phase of a free-falling atom, although of great physical importance, is technically much 
more difficult than the situation considered in our current paper. The reason is that in a 
Schwarzschild space-time, the vacuum fluctuations, namely, the environment the atom is 
coupled to, varies with the radial distance from the black hole. We will leave this to a future 
investigation. Nevertheless, as it has already been done in a recent experiment [26], we may 
manage to measure the decoherence factor of the off-diagonal elements of the reduced density 
matrix directly in order to get the geometric phase. So, we can obtain the geometric phase 
in a Schwarzschild space-time and that in a flat space-time respectively, and then compare 
the results in the two situations to reveal the effects of Hawking radiation. In summary, our 
analysis shows, at least in principle, that geometric phase might provide a potentially new 
way to detect the Hawking radiation, albeit difficulties of constructing specific experimental 
procedures that can actually be carried out. 
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